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Abstract 

The problem of velocity selection for reaction fronts has been intensively 
investigated, leading to the successful marginal stability approach for prop- 
agation into an unstable state. Because the front velocity is controlled by 
the leading edge which perforce has low density, it is interesting to study the 
role that finite particle number fluctuations have on this picture. Here, we 
use the well-known mapping of discrete Markov processes to stochastic dif- 
ferential equations and focus on the front velocity in the simple A + A —>■ A 
system. Our results are consistent with a recent (heuristic) proposal that 

VMS ~ V ~ j^. 
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I. INTRODUCTION 



There has been a great deal of interest in the problem of reaction front propagation in 
non-equilibrium systems. This issue arises in systems ranging from flames [l[ to bacterial 
colonies 0, from solidification patterns || to genetics M. Most of the theoretical work 
in this area involves solving deterministic reaction-diffusion equations. Here, we focus on 
effects that occur when one goes beyond this mean-field treatment and considers the effects 
of fluctuations. 

By now, it is clear that there are several possible mechanisms whereby the velocity of 
a deterministic react ion- diffusion front can be selected. For cases where we propagate into 
a linearly unstable state, the marginal stability criterion || suggests that the fastest stable 
front is the one that is observed for all physical initial conditions. For propagation into a 
metastable state, there is a unique front solution consistent with the boundary conditions 
and hence there is no selection to be done. In between, there is the case of a nonlinearly 
unstable state in which the exponentially localized front is chosen. These principles have 
been verified in many examples and in some cases can be rigorously derived ||. 

However, it is understood that deterministic equations are often only approximations 
to the actual non-equilibrium dynamics. This is particularly clear in the case of chemical 
reaction systems where the true dynamics is a continuous time Markov process which gives 
rise to a reaction-diffusion equation only in the limit of an infinite number of particles 
per site f7j. More generally, having a finite number of particles gives rise to fluctuations 
that may be important in the front propagation problem. It has been hypothesized in a 
variety of systems |8 Hl0fl that the leading effect of such fluctuations is to provide an effective 
cutoff on the reaction rate at very chemical concentrations. If this is the case, calculations 
by Brunet and Derrida JTT| predict that in the case of a system which (in the deterministic 
limit) exhibits (linear) marginal stability (MS) selection , the front velocity obeys the scaling 
Vms ~ v ~ h^lv' w ^ ere N is the (mean-field) number of particles per site in the stable state. 



Direct simulations of the underlying Markov processes have, in two cases to date |TT||T2]j , been 
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consistent with this predicted form, albeit with some uncertainty regarding the coefficient. 
Also, we note in passing that the cutoff idea is the simplest one which explains the recently 



discovered fact |I3| that one can have diffusive instabilities of a front in a chemical reaction 
system which do not show up in a reaction-diffusion treatment thereof. 

Our purpose here is to introduce a different approach for studying the role of these 
fluctuations in modifying the front velocity. There is a well-established machinery which 
transforms the master equation for Markov processes for chemical reaction systems to the 



solution of an associated stochastic differential equation. This was first proposed by Doi [|14 



and clarified in some seminal work of Peliti [15| . This framework has in fact been used for 



the study of critical phenomena associated with bulk transitions in reaction dynamics |T6 



but has not been applied to the issue of front propagation far from such a bulk transition. 
Here, we directly simulate the relevant stochastic equation; this requires the analytic solution 
of a (interesting in its own right) single-site problem, which then, via a split-step method, 
allows up to time-step the entire spatially-extended system. Our results to date verify the 
Brunet-Derrida scaling and in fact are even consistent with the coefficient obtained by the 
cutoff approach. 

The outline of this work is as follows. In section II, we review the mapping from the 
master equation to a Langevin equation with multiplicative noise. Next, we solve a variety of 
single-site problems as a prelude to introducing our simulation method. We then tackle the 
front problem numerically and compare our findings to the results obtained by augmenting 
the deterministic system with a cutoff. In order to accomplish this, the findings of Brunet 
and Derrida are extended to include the effects of finite resolution in space and time. At 
the end, we summarize the open issues that we hope to address in the future. 

II. DERIVATION OF THE STOCHASTIC EQUATION 

In this paper, we will study the following space-lattice reaction scheme: 

A^2A (I) 
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(2) 
(3) 



0^ A 



(4) 



(5) 



where e is the nearest neighbor sites of site i; a, /3, A, /i, r are rates of the corresponding 
reactions, i.e. probabilities of transition per unit time. This process is described by the 
master equation 



dP{{nA-t) 
dt 



= E 



dP({m};t) 



dt 



+ 



dP({ ni };t) 



dP({ ni };t) 



dt 

dP({n t };t) 



/3 + dt 

dP({n t };t) 



+ 



dt 



+ 



dt 



(6) 



which states that the probability P({rii};t) of having rii particles on sites % at some time t 
changes via each of the elementary processes: 

1. one particle splitting into two 

dP({m};t) 



dt 



= a[(rii - 1)P(. . . , rii - 1, . . . ; t) - UiP(. . . , m, . . . ; *)] (7) 



2. two particle reaction with one being annihilated 
dP({n t };t) 



dt 



p[(m + i)n;P(. . . , m + 1, . . . ; t) - mim - i)P(. ..,m,...;t)] (8) 



3. one particle annihilation 

dP({mY,t) 



dt 



X[(ni + 1)P(. . . , rii + 1, . . . ; t) - n,P(. . . , n h . . . ; f)] 



(9) 



4. particle birth from vacuum 
dP{{ni}-t) 



dt 



= r[P(...,n i -l,...;t)-P(...,n i ,. ..;*)] 



(10) 
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5. diffusion 

dP({n t };t) 



Of 



^ + 1)^0 ■■,n i -l,n e + l,...;t)- mP(. . . , n h . . . ; t)\ (11) 



In this section, we provide a self-contained derivation of the stochastic equation whose 
solution is directly related to the solution of this master equation. This is by now fairly 
standard, but we find it useful to include this derivation here both for completeness and for 
fixing various parameters in the final Langevin system. 

Following Doi [[L4}| , we introduce a vector in Fock space ({n,}) and raising and lowering 
operators af, hi with the properties: 

hi\ . . . ,n i} . . .) = 7ii\ . . . ,rii - 1, . . .) 
af\ . . . ,m, ...) = I ...,ni + 1, . . .) (12) 

and the commutation relation 

[a;,a+] = Sij (13) 

We choose an initial condition for the master equation (Q) to be a Poisson state, 

m 

P({n r };t = 0) = e- N ^Yl^ (14) 

i n i- 

where Na(0) = J2 n oi is the expected total number of particles. If we define the time 

i 

dependent vector 

l0(*)> = E p (W;*)l{"i» (is) 

{n{} 

the master equation can be written in the Schrodinger form 

J^(*)> = -H\m (is) 

where H = J2iHi and the latter is given by 

-/i^]a^(d e - hi) + a[l — hf]hfai — (3[1 — hf]hfhf — A[l — a^a,, + r[l — ojj"] (17) 



The formal solution of this equation is 

\(j)(t)) = e~ 6t \(f>(0)) = e - NA ^e~ At e^i & t noi \0) (18) 

To be able to calculate average values for observables, we need to introduce the projection 
state 

(i = (oin^ (19) 

i 

The external product of this with any state gives 1. Then any normal-ordered poly- 

nomial operator satisfies 

(\Q({at},{a l }) = (\Q({l},{a l }) (20) 
Using this equation we get for any observable 

(A(t)) = 5]il(W)P(K};() = (\%{&i})W)) (21) 

{rii} 

where A({di}) is what we obtain by using the commutation relation to normal order A and 
thereafter setting af to 1. 

In order to write a path integral representation for the time evolution operator, we 



introduce a set of coherent states ( see |17]] for more strict treatment) 



\\{zi}) = <£'«*t- Ml \0) (22) 

where z% is complex eigenvalue of cij. In the case of real positive {z{\ this states are Poissonian 
states with (rij) = Zi. By inserting the completeness relation, 

1 = (IX / — e-l^+^HzX^I (23) 

where d 2 Zi = d(Rezi)d(Imzi) , into expression (pH]) we get 

(|i({a,})e-^|0(O)) = (|I({a i })(l - HAt) NAt \\z^} = 

j2 ( n ) 

(\Mm)([[ I ^- e -i^ ) i 2 +^ ) +< w )pw)(zW|i x 

^ d 2 z\ 

j=l j=l i 
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d 2 ^ 

TV 

N N , J2 Jj) 



:n(i-^A0)ii^> = (oi(nn / ^h^^d^io) m 



where At = t/N, \<f>(0)) = \\z^) and 

TV 

S = E E( z^)At + \z^\ 2 - zfht^ - 4 N) + 4 0) ) = (25) 

j=l i 

E EM- — ^7 — + + 1; ^ (J))) (26) 

j=0 i 

Here = zf ] + 1 and Hi({z* {j+1) }, {z®}) is the same function of zf j+1) , z* U) as 

Hi({af }, {a,i}) of af , a,i. In the continuous time limit we get 

(Aft v - ;n^^^({^fa)})e^ [{z -' (f)} ' {2 ' (f)} ^ (97 , 

[ {0)) JlliVziVzie- 3 ^}^)}^ 1 ' 

with 

S[{zi(t)}, {*(*)}; t ] = E jf° dt (^(*)[^ ~ AiV 2 ]^(t) - a(l + z i (t))z i (t)z i (t)+ 

(3(1 + Zi(t))zi(t)z?(t) + \zi{t)zi{t) - rzi{t)) (28) 



where V 2 is the lattice Laplacian Zi(t) = J2 e { z e{t) — Zi(t)). Now, we linearize the action 
using the Stratonovich transformation 



e m«,-H>^ j Ahe -i^v*=^m^ (29) 

and integrate out the z variables 

J i t=o 

azi(t)dt + Pz?(t)dt + Xzi(t)dt - rdt - ^2(aZi{t) - (3z^(t))rji(t)Vdt) (30) 



In this expression, there are 5-functions at every time; this means that only Zi(t) which 
satisfy to the Langevin equation 



dzi(t) = fiV 2 Zi(t)dt +{a- X)Zi(t)dt - (3zf{t)dt + rdt + y/2(aZi(t) - /3z?(t))dWi(t) (31) 

(where Wi(t) is a Wiener process) contribute to the path integral. In other words, the 
variables Zi(t) remain on the trajectories described by equation fl3~T|). Note that this equation 
must be considered as an Ito stochastic differential equation, since we can see from the form 



of the action ( pq) that the updating the variables Z{ to time-step j+1 only requires knowledge 
of the variables at time-step j. Also, we note that for A > and for small enough (positive) 
r, if the initial conditions specify < ^(0) < a/ (3, this will remain true for all subsequent 
time. Thus, equation fl31| ) describes the temporal evolution of the system as sequence of 
Poissonian states [TE| . 



For further analysis, we rescale (|3lj) with z = ua//3, t — > t/a, A = A/a, f = rj3/a 2 , 
and fi = n/a. If we furthermore let N = a/ (3 be the mean-field number of particles in the 
presence only the first two processes (no spontaneous decay or spontaneous creation), we 
obtain 



du-i = jT\7 2 Uidt + (1 — \)uidt — u 2 dt + fdt + y J^V Ui ~ u fdWi (32) 
with initial conditions < -Uj(O) < 1. 



III. EXACT SOLUTIONS OF SOME LOCAL LANGEVIN EQUATIONS 

In the absence of process (f|), i.e. at r = 0, equation (^) has an absorbing state u = 0. 
In the vicinity of this point, equation (|32D cannot be treated by merely setting dt to a finite 
time-step. Such a scheme would often give rise to a negative u, due to the (very-large) 
noise term. One ad- hoc way to circumvent this difficulty was given by Dickman |19| , who 
proposed to re-introduce discreteness into the state-space in the vicinity of the absorbing 
state. Although this approach appears to work (it seems to lead to the correct critical 
behavior near the bulk transition of this class of models), it seems to be a step backward; 
after all, the original process was discrete and the whole purpose of using the Langevin 
formalism is to provide a (hopefully more analytically tractable) continuum description. 
But, one must then come up with a different scheme for updating the stochastic variables. 

Our approach is to solve exactly the stochastic part of the evolution equation and embed 
this via the split-step method in a complete update scheme for a finite time-step At. We 
will discuss the details of this scheme in the next section. Here, we provide an analytic 
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solution for several (local) Langevin equations, as these results will be needed later. Also, 
this solution set is of interest on its own. There is some limited consideration of equations 
of this sort in the literature ||20|| , but as far as we can determine, these explicit solutions for 



the case of physical no-flux boundary conditions at the absorbing state has not previously 
appeared. 

So, we consider Langevin equations with just the noise term. Let us start with the 
simplest example, 

du = VZudW (33) 

The probability density P(u, t) satisfies the associated Fokker-Planck equation 

dP(u,t) d 2 „. , 

at = u ' *' < 34) 

with initial condition P(u, t) = 5{u — u ). We want our solution to be equal to zero at 

t — tin 

u < and have no flux leaking out of this point; this will guarantee that the total probability 
remains a constant, which we will choose to be unity. 

To solve this equation, we define ip = uP and Laplace transform in time to obtain 

d 2 ip 

Here, ip(s) is the transform of ip. If we let y = 2y/u, this can be written as 

(36) 

The homogeneous part of this equation can be recognized as a variant of Bessel's equation. 
This allows us to write down a provisional solution in the form 

$( 8 ) = ^K x Whl>)hWh<) (37) 

where I\ and K\ are modified Bessel functions and y > (y < ) is the larger (smaller) of y and 
yo- Returning to the original variables, 



P{s) = 2^if 1 (2 v /i^)/ 1 (2 v /^) (38) 

\/U 



This solution does not, of course, vanish for u < and hence we must modify it by multiply- 
ing by 9(u). This does not change the fact that it solves the equation away from u = but 
it does introduce a discontinuity of size 2 ^JuqK\ {2^/suq) . If we look at the original equation, 
we see that this leads to a 5 function via u8'(u) = —S(u). This must be compensated by 
adding an explicit S function piece to the solution. The final result is 



P(s) = 2^ J fir 1 (2 v /i^)/ 1 (2 v /5^) + 2J—K 1 (2y/au^6(u) (39) 
\/u V s 



We can do the inverse transform by the usual contour integral approach. The details are 
particularly unilluminating, so we merely quote the final result 

P(M ) = I &-(«+«o)A /i( V^o ) + e -«o/ i(5(M) (40) 

t V UL t 

One can check explicitly that this solves the equation and also that P remains normalized 
for all times. The 8 function piece represents accumulation at the absorbing state; as t gets 
large, all the probability end up there. The regular part of P(u, t) is presented on the Figure 
0. We see that as ratio uo/t becomes smaller, the distribution gradually shifts towards zero 
and differs from the Gaussian expected at very short times. 

For completeness, we write down the solution of Langevin equations with additional 
terms. If we take the system, 



du = rdt + V2udW, (41) 

the probability density is 



p( M ) = ' t > e -(u+u )/t (42) 

For this case, spontaneous birth from the vacuum prevents the system from falling irre- 
versibly to the state u = 0. Instead, there is an integrable power-law singularity near u = 
which becomes a 5 function in the r — > limit; this is shown in Figure |[ 
For the system 



du = audt + V2udW (43) 
10 



the probability density is 

P(u,t) = e-^-^(tt) + \jg{t)^h (2yjg{t)uu ) e~ a ^ (44) 

where g(t) = / ■ 111 this case, the drift toward infinity gives rise to a finite total 
probability (for long times) of falling into the absorbing state. One can also work out the 
case of both finite r and finite a. 
For the system 



du = ^2{u - u 2 )dW (45) 
we can derive a series representation for the probability density, 



+°° 9m -I- 1 

P(u,t)= £ e -™(™+i)* 

m=l 



2m(m + 1) \ u(l 



m \ 

U 



-OO 



2m + 1 



*(«) [1 - «o - g (-l) W 2^pi) V"°(l - «o)e-^ +1 ^P^ t (2 Mo - 1)] 



+oo 



2m + 1 



5(1 - «)[«o + £ o 7 M. N 1 - «o)e- m(m+1) ^(2 Mo - 1)] (46) 



rn= 



where P^x) is Legendre polynomial. Successive terms in this sum decay rapidly because of 
the fast exponent, and the sum can be computed numerically to high accuracy. Note that 
there are absorbing states at both u = and u — 1. If the initial state starts close to one of 
these, the probability is almost the same as (fTO|); if the initial state is in between, then for 
short times the system is almost Gaussian. 



IV. FRONT PROPAGATION 

A. Numerical method 

We are interested in numerically solving equation (|3l"D , for the particular case of r = A = 
0; this is case which reduces in the deterministic case to the well-studied Fisher equation || 
For this purpose, we define a function 

F U0 ,At(u) = f U _P(u,At)du (47) 

J 
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where P(u, At) is the analytical solution of a single-site Langevin equation such as fl32|). This 
function has values ranging from to 1. If y is random variable homogeneously distributed 
on [0,1], then u = F~ oAt (y) is distributed according to corresponding truncated Langevin 
equation at time At. The remaining parts of complete Langevin equation are deterministic 
and for those terms we can update u via a simple Euler scheme. We then can combine this 
two steps together; we first compute the change in u due to fluctuations and then the change 
of u due to the deterministic part (using new value of u). Thus 

u(t + At) = F- { l )At (y) + V{F- ( l )At (y)}At (48) 

where T>{u] denotes the terms remaining after consideration of the noise term. 

It is important to note that this scheme never allows the field to go below zero, but 
it does allow a variable to be stuck at zero until it is "lifted" by the diffusive interaction. 
This is an absolutely necessary aspect of simulating processes with an absorbing state. 
Approximations which do not allow for getting "stuck" , such as the system-size expansion 
method of Van-Kampen ]7J (where the noise correlation is taken to be related to the solution 
of the deterministic limit of the equation) get this wrong and hence cannot get the correct 
front velocity. This explains why the simulation results of [EHl] do not at all exhibit the 



anomalous N dependence expected via the Brunet-Derrida cutoff argument. As we will see, 
our approach is much more successful. 



B. Marginal stability criterion for a discretized Fisher equation 

As N gets large, our results should approach those of the deterministic system. Since 
this problem corresponds to propagation into an unstable region, the velocity should be 
given by the marginal stability approach. As is well known, this predicts a velocity equal to 
1\f~D 1 in the continuum (in time and space) limit. Here, we extend this result to a discrete 
lattice and a finite time update scheme, so as to be able to directly compare our simulation 
data with the theoretical expectation. 
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After linearization, the deterministic part of the discretized equation takes the form 

H(n^i-2u? + u^ 1 )+u^ (49) 



At 

We want to compare this equation to the usual Fisher equation with diffusion coefficient D. 
This means that /i = D/h 2 ; we will consider the case D = 1. We assume that the front 
moves with constant velocity c and therefore the variables show a stroboscopic picture 
of this motion at times j At on the lattice sites i. If we move with the speed of the front we 
will see that its shape exponentially decays as e q<y%h ~ c ^ At \ Substitution of this expression to 
(IH) gives the dependence of c on the decay rate q 

e -qcAt _ i 9 

— = _(cosh^-l) + l (50) 

The standard marginal stability argument predicts that we can determine the decay rate 
and asymptotic speed of the front (for a sufficiently localized initial state) by solving (|50"1) 
as well its derivative with respect to q 

ce -gcAt = 2 . nhqh ^ 

1 1 



Simulations directly confirm this formula as well as the Brunet and Derrida ||TT[ result 
(actually derived earlier by Bramson that — c(t) ~ 1/t (see Figures |3|). 

As already mentioned, it has been conjectured that the leading effect of the fluctuations 
is the imposition of an effective cutoff of order 1/N in the deterministic equation. To check 
this, we need to extend the Brunet and Derrida result to our discretized equation. The 
basic idea is that there must be a small imaginary part of the decay rate so at to satisfy 



the continuity conditions at the cutoff point; this is discussed in detail in pT| . This leads 
directly to 

where go is solution from the marginal stability criterion and A is some constant. Since c'(q) 
is zero at the marginal stability point, we can find the change in velocity 5c by considering 
the second derivative of the function c(q) given by (H). We thus get 
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Again, simulations confirm this formula (see Figure |4j). 



C. Results 

We now present the results of our simulation. We chose to make one further simplifi- 
cation. We use the pure square root noise term instead of the precisely correct term given 
in equation ([32]). We do this for computational ease, inasmuch as the expression derived 
for this case is much simpler than that of equation (p5|). Since it is only the effect of the 
noise near the u = absorbing state which is crucial for altering the selected velocity, this 
simplification should not be essential. Once we have done this, the resulting equation has 
the nice feature that the coefficient 1/y/N in front of the noise term can be removed by the 
re-scaling u = uN . This means that we can simulate equation ( |32D using a fixed probability 
table (with the same time step) for any N. 

To actually evaluate the probability table F U0) & t (u), we chose 512 equidistant values of 
uq in the interval from to 30. For each u , the interval of values for u where F„ 0i At=o.oi(w) is 
non-trivial was divided into 1024 equidistant points. The new value of u was then determined 
by linear interpolation of the data from the table. For uq > 30, new values of u were 
determined using a standard algorithm for the Gaussian distribution , since this distribution 



is the asymptotic limit of equation (4"0), when u 3> At, u ^> At. The difference for 
this distribution and exact solution is small for uq > 30. Finally, the computation of the 
stochastic term was turned off for uq > 10~ 3 iV. This should not affect the speed which, we 
have already argued, is only sensitive to what happens near u = 0; this insensitivity was 
also checked directly by running some simulations in which the stochastic term is included 
for all values of u. 

All of our simulations were run up to the time t = 10 6 . Four values of the velocity, 
corresponding to time intervals of approximately 2.xl0 5 were obtained so as to get an average 
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and some error bar. In Figure |5|. we show data in the form of c — Cqo, where Cqo is calculated 
from equation (|5T|), versus In 2 N. Also plotted for comparison is the function ll/ln 2 iV. 
Note that over many orders of magnitude of N, the dependence derived by Brunet and 
Derrida provide a very good fit to the data. 

To get a more accurate indication of the data for large N, we present in Figure |6| a series 
of three runs, for differing values of the spatial lattice spacing h. Under the hypothesis that 
the stochastic system should be precisely the same as the deterministic system with the 
cutoff added, the expected limiting values at infinite N are shown as triangles on the axis. 
It is clearly impossible to definitively conclude that the curves are approaching these values. 
On the other hand, simple extrapolations come very close and we believe that it is more 
likely than not that this hypothesis is true. This is opposite to what was conjectured based 
on simulations of a discrete Markov process, where the velocity seems to scale albeit with a 
different coefficient. Given the incredibly slow convergence of this velocity at large N, we 
are pessimistic as to whether any purely simulational strategy would provide a definitive 
answer to this question. This therefore offers a crucial issue for future theoretical analysis 
to investigate. 



V. SUMMARY 

In this work, we have shown how to use the field-theoretic mapping of discrete Markov 
processes to stochastic equations for continuous density variable to address the role of finite- 
particle number fluctuations on the velocity of reaction fronts. Specifically, we studied a 
model which leads to the well-known Fisher equation in the N — > oo limit, where N is the 
average number of particles per site in the stable state. Our goal is to understand how the 
usual marginal stability criterion becomes modified by these stochastic effects. 

It is clear that having finite N lowers the velocity at which a front (corresponding to the 
invasion of the unstable state by the stable one) will propagate. One attractive hypothesis 
is that the leading effect of the fluctuations is to to introduce an effect cutoff into the 
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deterministic equation; this idea arose independently in model of biological evolution || and 
in mean field approaches for diffusion-limited-aggregation (DLA) Brunet and Derrida 
have shown that if this is the case, one should expect vms — v — j^y, where C = tt 2 for 
the case of continuous time and space. We have extended the calculation of C to the finite 
lattice size, finite time-step system and compared this prediction with direct simulations of 
the relevant stochastic equation. Our results verify the form of the scaling and suggest that 
the coefficient may be correct as well. 

One issue that is left unaddressed by our work to date concerns the effects of higher 
spatial dimensionality. It is likely, although unproven, that the velocity change will be 
smaller, as the fluctuations get averaged over the transverse directions. This seems to be 



the explanation for the findings of Riordan et al p3| that the reaction front looks mean-field 
like even for small N, in three and four dimensions. We hope to report on this issue in the 
future. 

Finally, we point out yet again that there is no analytic treatment available for the 
velocity selection problem in the stochastic equation. Obvious expansion methods such as 
the system-size approach cannot work, as they neglect the essential role of the fluctuations 
to push the system back into the absorbing state at small density. We need to find a more 
powerful approach! 

We thank D. Kessler for many useful discussions. We also acknowledge the support of 
the US NSF under grant DMR98-5735. 



16 



REFERENCES 

[1] A. I. Kolmogorov, I. Petrovsky and N. Piscounov, Moscow Univ. Bull. Math. 1, 1 (1937). 

[2] I. Golding, Y. Kozlovsky and I. Cohen and E. Ben- Jacob, em Physica A, in press (1998) 
and references therein. 

[3] For a review, see D. A. Kessler, J. Koplik and H. Levine, Adv. Phys. 37, 255 (1988). 

[4] R. A. Fisher, Annual Eugenics, 7, 255 (1937). 

[5] E. Ben- Jacob, H. Brand, G. Dee, L. Kramer and J.S. Langer, Physica D, 14, 348 
(1985); W. Van Saarloos, Phys. Rev A 39, 6367 (1989); G. C. Paquette, L.-Y. Chen, N. 
Goldenfeld and Y. Oono, Phys. Rev. Lett. 72, 76 (1994). 

[6] D. G. Aronson and H. F. Weinberger in Lecture Notes in Mathematics; Partial Differ- 
ential Equations and Related Topics , A. Dold and B. Eckmann eds., Springer- Verlag 
(Berlin, 1975). 

[7] N. G. van Kampen, Stochastic Processes in Physics and Chemistry (North-Holland, 
Amsterdam, 1981) 

[8] E. Brener, H. Levine and Y. Tu, Phys. Rev. Lett, 66, 1978 (1990). 

[9] T.B. Kepler and A.S.Perelson, Proc. Natl. Acad. Set. USA, 92(18), 8219 (1995). 
[10] L. Tsimring, H. Levine, and D. A. Kessler, Phys. Rev. Lett. 76, 4440 (1996). 
[11] E. Brunet, B. Derrida, Phys. Rev. Em, 2597 (1997). 
[12] D. A. Kessler, Z. Ner, L. M. Sander, Phys. Rev. £58, 107 (1998). 
[13] D. A. Kessler and H. Levine, Nature 394, 556 (1998). 
[14] M. Doi, J. Phys. A 9, 1479 (1976). 
[15] L. Peliti, J. de Physique 46, 1469 (1985). 



17 



[16] J. L. Cardy, U. C. Tauber, J. Stat. Phys. 90, 1, (1998) and references therein. 
[17] M. Cafaloni, E. Onofri, Nuclear Phys. B 151, 118 (1979). 
[18] G. W. Gardiner, Handbook of Stochastic Methods (Springer, Berlin, 1985). 
[19] R. Dickman, Phys. Rev. £50, 4404 (1994). 

[20] W. Feller,v4rm. Math. 54, 173 (1951); H.-J. Buttler and J. Waldvogel, Mathematical 
Finance 6, 53 (1996); J. Pitman and M. Yor, Z. Wahrschein. verw. Gebiete 59, 452 
(1982). 

[21] A. Lemarchand, A. Lesne, and M. Mareschal, Phys. Rev. E 51, 4457 (1995). 
[22] M. Bramson, Mem. Am. Math. Soc. 44, No. 285 (1983). 

[23] J. Riordan, C. R. Doering and D. Ben-avraham, Phys. Rev. Lett. 75, 656 (1995). 



18 



FIGURES 




19 




20 



1.66 i 1 1 1 

1.64 

1.62 - • • * " . . ..• ••.* • 

a . -„. • 

l • • 

8 • *• 

3, . . • • 

1.6 ' ' ' - 

1.58 

1 .56 1 1 1 

5000 10000 15000 

t 

FIG. 3. Front propagation for the discretized Fisher equation. Parameters were At = 0.01, 
h = 1; for this choice, equation (|49|) gives Cqo — 2.054115884. The speed is determined by registering 
every time some lattice point comes within some small interval around .5 and then merely using 
the ratio of the number per spaces moved over the time elapsed; this method leads to some scatter, 
as seen in the graph. Aside from showing agreement with the calculated Cqo, the graph also shows 
how the system approaches the eventual 1/t in asymptotic behavior. 
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FIG. 4. Dependence of the front speed on the deterministic cutoff. The curves represent simu- 
lation results calculated with the parameters h = 1 and At = 0.01; 0.015; 0.02 from top to bottom 
respectively. The cutoff was implemented via setting the field to zero at any time and place where 
it was below 1/N. Points on y-axis represent theoretical values when N — > oo as given by equation 
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FIG. 5. Front speed for the Langevin equation with parameters At = 0.01 and h = 1. 
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FIG. 6. Front speed for the Langevin equation with parameters At = 0.01 and h = 1, h = 1/2, 
/i = 1/3 correspondingly from top to bottom. Triangles on the y-axis are the corresponding 
asymptotic values for the deterministic cutoff. 
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